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Abstract 

The identity of the quark-squark-gluino Yukawa couphng with the corresponding 
quark-quark-gluon QCD couphng in supersymmetric theories can be examined experi- 
mentally at the Large Hadron Collider (LHC). Extending earlier investigations of like- 
sign di-lepton final states, we include jets in the analysis of the minimal supersymmetric 
standard model, adding squark-gluino and gluino-pair production to squark-pair pro- 
duction. Moreover we expand the method towards model-independent analyses which 
cover more general scenarios. In all cases, squark decays to light charginos and neu- 
tralinos persist to play a dominant role. 



1 Introduction 



A characteristic consequence of supersymmetry is the identity of the quark-squark-gluino 
Yukawa couphng with the quark-quark-gluon couphng in SUSY-QCD [1]: 

9s{qqg) = 9s{qqg) (1) 

which persists to hold after the supersymmetry is broken softly. Similar identities are pre- 
dicted in other sectors of the supersymmetric extension of the Standard Model. 

Pair production of sleptons and charginos/neutralinos in e+e^ collisions has been inves- 
tigated to test the identity of Yukawa with gauge couplings in the non-coloured electroweak 
sector, cf. Refs. [2-4]. The equality of electroweak Yukawa and gauge couplings is expected 
to be studied at the per-cent to per-mil level by these experiments. In a recent report, 
Ref . [5] , the measurement of like-sign di-leptons has been examined to study the squark-pair 
production processes qq —>■ qq. Since these processes are mediated solely by the exchange of 
gluinos, they provide, in principle, the most appealing measurement of the Yukawa coupling 
in supersymmetric QCD. All branching ratios in the decay cascade <? — > Xi^ ^ have been 
assumed to be known in evaluating the qq pair cross section which involves the Yukawa 
coupling. 

In this note we extend the analysis in two directions. First we will include the measure- 
ment of jets in like-sign di-lepton events, motivated by the large production cross sections for 
gluinos [6] which significantly enlarge the like-sign di-lepton ensemble through instantaneous 
decays to real or virtual squarks [see also Ref. [7]]. In a second step we will investigate 
experimental scenarios in which only the properties of light gaugino-type charginos and neu- 
tralinos have been pre-determined either at LHC or at an e~^e~ linear collider [8]. This 
method, within a wide variety of theories, does not rely, in particular, on the measurement 
of the decay branching ratios of heavy squarks. Instead of studying the absolute size of the 
cross sections, we will analyse ratios of like-sign di-lepton cross sections and the evolution 
of these cross sections with the transverse momenta of additional jets. Including jets with 
varying transverse momenta in the analysis changes the relative weight of the subprocesses 
with different dependence on the Yukawa coupling. Therefore, relating the cross sections 
for different final state configurations, the overall normalisation can be eliminated from the 
analysis and the Yukawa coupling can be determined without a priori knowledge of the 
branching ratios for q decays. 



2 Squark and gluino production at the LHC 

Most suited for this purpose, in the area of parameter space in which gluinos are heavier 
than squarks, are the basic subprocesses 

squark-pair production '■ qq—>-qq{g) (2) 
Compton process : gq qg (g) with g ^ qq (3) 

cf. Fig. [H Gluinos decay instantly to squarks under this kinematical condition, {g) denotes 
additional gluon jets emitted in the production process. The analysis of this kinematical 
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Figure 1: Characteristic examples of Feynman diagrams for partonic squark and gluino 
production in hadron collisions. Dots indicate the gauge coupling gs, while squares represent 
for the Yukawa coupling Qg. Gluon radiation in the initial state, indicated by the dashed 
lines, can lead to extra hard jets. 

area, admittedly, is easier than the reverse mass configuration but sufficient to illustrate the 
basic points. Other processes of secondary importance, like gluino-gluino pair production 
etc., have been included in the final analysis properly. As evident from Fig. [H the size of 
the cross section for squark-pair production is set by the fourth power of the strong Yukawa 
coupling, g^, while the cross section of the Compton process is determined by the product of 
the second power of the gauge and Yukawa couplings, g^gs- Since the relative weight of the 
subprocesses changes for different jet ffnal-state conffgurations, the Yukawa coupling can be 
measured by comparing the cross sections in different parts of the phase space. 

To focus on the basic processes Eqs. ([2]), ([3]) and to reduce the backgrounds from other 
SUSY processes, squark-squark pairs of the same charge should be selected. This can be 
achieved by studying like-sign di-lepton final states [5] in gaugino-type chargino and neu- 
tralino decay chains of two L- and R-type squarks like 

Ul ^ d + xt jet + £++^ 
WL,R u + xl ^ jet + £+/r +^ ^ ' 

dh u + Xi jet + r + ^ 

4r ^ d + xl jet + r/£++^ 

with ^ denoting the missing energy. SU(2) isospin invariance guarantees that the partial 
widths for ml and c/l decays to charginos are equal to very good approximation in this 
configuration. Assuming CP-invariance, the subsequent leptonic decays of Xi ^^^^ related 
to xf, and the g£ decays to the decays. For pure SU(2) singlet gaugino LSP states Xi; 
R-type squarks will decay into the channel qn + Xi^ iiot generating charged leptons in the 
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final state. 

Extended mixing in the third-generation stop and sbottom sector will significantly in- 
crease the complexity of the analysis [though the production of stops may prove particularly 
useful in testing the Majorana character of gluinos [9]]. However, scalar top and bottom de- 
cays generate b quarks in the final states which, by tagging, can be used to strongly reduce 
the influence of these sparticles on the analysis. 

The branching ratios for gluino decays to squarks plus quarks depend only on the gluino, 
squark and quark masses so that they can be predicted unambiguously in the Compton 
process. To leading order the widths are given by the Yukawa couplings which however 
cancel from the branching ratios. 



3 Di-lepton plus jet analyses in the MSSM 

In the first step we have elaborated on the details of the scheme outlined above within the 
SUSY scenario defined in Ref. [5]. Apart from the slightly increased gluino mass to improve 
on jet tagging, the scenario corresponds to the reference point SPSla [10], describing a spe- 
cific light-mass scenario within the minimal supersymmetric standard model MSSmIiI Typical 
parameters, masses and branching ratios relevant to the present analysis, are collected in 
Table [D 

Similarly to Ref. [5], we will assume in this section that all the parameters, masses and 
branching ratios, have been pre-determined in squark cascades at LHC and/or analyses 
at an e"'"e~ linear collider. It should be noted, however, that within the group-theoretical 
frame of the minimal supersymmetric extension of the Standard Model [MSSM], the squark 
branching ratios can be predicted if the structure of the chargino/neutralino sector has been 
explored at a low-energy lepton collider. Within this framework the measurement of the 
squark-decay branching ratios can be substituted by the completeness assumption of the 
MSSM. The assumption can be cross-checked internally at LHC by searching for non-MSSM 
squark decays. 

About two thirds of the time, the light-flavor L-squarks decay into the charginos x^, 
while about one third of the decays go into the neutralino X2, following from the fact that 
the L-squarks have small hypercharges, Yg^ = 1/6, and thus predominantly couple to winos. 
Furthermore, since the mixing between the bino and wino components of the neutralinos 
is small, as well as the mixing between L- and R-squark states negligible, the branching 
fractions of the up- and down-squarks are almost equal. This is in particular the case in 
our reference scenario, cf. Table [H and a consequence of SU(2) and CP invariance. [Note 
that in general CP invariance can be broken in the MSSM by complex parameters, but in 
this analysis all supersymmetry and supersymmetry breaking parameters are assumed to be 
real.] 

Due to the relatively large value of tan/3 = 10 in this scenario, the charginos al- 
most exclusively decay to the scalar tau fi, which subsequently decays to taus. To trace 

-'^The model in Ref. [5] is defined strictly at the weak scale and is not intended to accurately reflect either 
GUT scale physics or cosmological parameters which most likely are affected by other mechanisms too. 
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a) Weak scale soft supersymmetry breaking parameters: 
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c) Relevant branching ratios: 

BR[ui^^ dxt]=Q5% BR[x^ f^z/r] = 95% BR[g^qqj^]= 7.1%, q^b,t 

BR[dL ^ M xr] = 61% BR[x^ ^ W*^Xi] = 5% BR[^ ^ g gn] = 9.6%, g 7^ 6, t 

BR[ml ^ u xl] = 32% BR[^ -^bbi] = 23% 

BR[(iL ^ X2] = 31% BR[g t ii] = 10% 

Table 1: Reference scenario used in this analysis. It coincides with the MSSM Snowmass 
point SPSla [10], except for the shghtly larger gluino mass [5]. The low-energy spectrum 
has been calculated by means of SoftSUSY 1.9.1 [11]. 

the charge of the squarks exphcitly, wc restrict ourselves to the leptonic tau decays. The 
phenomenological analysis we will therefore be designed such that the decay chains 

65% , 100% , ,-, 35% 

Hl > dxt ^dr^UrXi ^df.^ + ^ 

61% _ 100% _ _ „ 35% _ „ r 

(Zl uxi ^ ur UtXi ^"^-^ V^ — ^tIA (6) 

are singled out. The numbers above the arrows denote the branching fractions. 

The LHC processes which have been studied in this report are hke-sign di-lepton final 
states plus a number of jets with variable transverse momenta: 

pp + 2 jets + $ (7) 

pp e^e^^ + 3 jets + ^ (8) 
pp e^e^ + 4 jets + ^ (9) 

and possibly additional jets from gluon radiation, that however tend to have smaller trans- 
verse energies than the jets emitted in squark and gluino decays. 
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The first process, Eq. (^^, concentrates on direct squark-pair production, the classical 
reaction [5] for measurements of the strong Yukawa couphng with the basic cross section 
scahng as g^. The transverse momentum of the jets is assumed beyond 200 GeV. 

In the second process, Eq. ([8]), we demand two high transverse momentum jets with 
Et > 200 GeV, plus an additional third jet with variable transverse momentum. For low 
transverse momenta of the third jet, direct squark-pair production is projected out, moving 
the dominant weight to the Compton process however for increasing transverse momentum. 
In this way, depending on the transverse momentum, different powers of the Yukawa coupling, 

~ Ci(pt)5's + C2(pt)5's^s; become relevant. 

The evaluation of the third process, Eq. ([9]), parallels the previous analysis. It does not 
provide new methodological insight but increases statistics when included. 

To analyse the prospects for extracting the SUSY-QCD Yukawa coupling in a reahstic 
experimental environment, we have generated samples for the squark and gluino production 
process as well as SM backgrounds with Pythia 6.408 [12], including gluon and quark 
radiation through virtuality-ordered parton showers, string hadronisation and the underlying 
event according to 'Tune A' [13]. The required modifications for varying the SUSY-QCD 
Yukawa coupling independently of the gauge coupling were implemented in a private versioij^, 
according to the lowest-order cross sections given in [6]. The main SM backgrounds originate 
from tt production as well as WWjj production. The WWjj process was generated with 
MadGraph [14], without QCD radiation and fragmentation. 

The cross-sections for squark, gluino and top production were normalised by i^"-factors 
to the values obtained with next-to-leading order QCD corrections [6, 15], while for the 
W^W^jj background only leading order results are available. In analysing distributions, 
K-factors as calculated for the total cross-sections have been included; this approximation is 
based on conclusions of Ref. [6] where the primary distributions have been shown not to be 
distorted by more than 10% so that a uniform rescaling is justified. A possible complication 
is created by the fact that the gluino is quite wide in our scenario, with Vg ~ 35 GeV at 
leading order, generating doubts whether the narrow- width limit is justified. For example, 
the convolution of the corresponding Breit-Wigner function(s) with the steeply falling parton 
distributions greatly distorts the final distribution - to the point of even creating a spurious 
peaking towards low masses. We note that this issue is generic for gluon-induced processes 
involving wide resonances at the LHC, and although we do not aim to solve it here, the 
reader should be aware of this problem. 

Finally, the detector response was parametrised by simple acceptance cuts: a toy calorime- 
ter spanning the pseudorapidity range jr^l < 5 with a resolution ofO.lxO.linr^x^ space 
was assumed. Jets were identified with a simple UAl-like cone jet algorithm (PYCELL), with 
a cone size of AR = 0.4. Muons and electrons were reconstructed inside \7]\ < 2.5 and an 
isolation criterion was further imposed, requiring both less than 10 GeV of additional energy 
deposited in a cone of size AR = 0.2 around the lepton and also no reconstructed jets with 
Ptj > 25 GeV closer than AR = 0.4 around the lepton. These criteria duplicate the default 

^Publicly available from Pythia version 6.4.19. 
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Process 


Total cross-section 


Cross-section with BRs and cuts 


Efficiency 




2.1 pb 


6.1 fb 


0.29% 




1.4 pb 


3.1 fb 


0.23% 




7.0 pb 


7.6 fb 


0.11% 


(^[99] 


3.2 pb 


1.4 fb 


0.04% 


a[SM] 


800 pb 


<0.6 fb 





Table 2: Signal and background cross-sections after applying the cuts mentioned in the 
text. 

settings of the ATLFAST simulation package [16]. 

All the signal processes, Eqs. ([7]) to ([9]), are characterised by the similar features of two 
leptons, missing energy and at least two hard jets. Therefore the following preselection cuts 
have been used in the event generation: 

• At least ^ > 100 GeV of missing energy. 

• At least 2 jets with pxj > 100 GeV. 

• Exactly two isolated same-sign leptons i = e, fi with p^/ > 7 GeV. 

To reduce the SM background, the final analysis in addition employs the soft selection cuts 

• At least ^ > 300 GeV of missing energy. 

• At least 2 jets with pxj > 200 GeV. 

• Veto on identified bottom-quark jets, using a tagging efficiency of e = 90% with a 
mistagging rate D = 25% [5, 16]. 

These cuts have been optimised to minimise the remaining SM background, however at the 
cost of being relatively expensive on the signal statistics as well. This choice is justified 
by the goal to interpret the measured signal rate and distributions as a combination of u, d 
squark and g gluino production processes, so that it is mandatory to keep the contamination 
from background processes as low as possible. 

The veto on bottom quarks is effective against the ti background, but it also helps to 
reject the unwanted decay modes of the gluinos into scalar bottoms and scalar tops. Since 
the stop and sbottom decay chains can deviate significantly from those of the light-fiavour 
squarks, it is advantageous to remove them from the analysis. As mentioned above, the 
relevant branching ratios of the gluinos into different squark flavours can be determined 
unambiguously from the gluino, squark and quark masses. 

After application of the cuts, the signal and background cross-sections listed in Tab. [2] 
are obtained. Thus, with a total luminosity of 300 fb~^ about 2700 squark pair events and 
1750 events from the Compton process are expected. 
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a[fb] 




Figure 2: Results of simulations for the number of identifiable jets, for the different processes 
that lead to the general signatures in Eqs. ([7]) to ([9]) . For the first two jets, a threshold of 
PTjet > 200 GeV has been used, while pxjet > 50 GeV for the other jets. 



After reducing the SM background, the three high-energy processes pp — > gLQ'L, QhSi, 99 
should be distinguished in order to extract the SUSY-QCD Yukawa couphng ^s- As outhned 
above, in a first step this can be achieved, statistically, by analysing the number of hard jets 
in the final state, see Fig. [2l With a transverse momentum threshold of PTjct > 50 GeV for 
the jets beyond the two hardest jets (with Pxjet > 200 GeV), the sample of events with only 
two jets is dominated by squark pair production. On the other hand, events with three or 
more jets largely originate from the Compton process pp — > qj^g. For A'jets > 4, gluino pair 
production also contributes to the di-lepton + jets signal, though the rate for this process is 
strongly suppressed due to the b veto. 

From the measured rates for events with 2, 3 and 4 jets, and using the pre-determined 
branching ratios for squarks and gluinos, the coupling can be determined. A similar 
analysis has already been performed in Ref. [5], but only using the 2-jet portion. Including 
the higher statistics of the 3- and 4-jet events, we conclude that a much more precise mea- 
surement of could be possible. From a fit to the simulated data, the statistical error is 
estimated to be only 0.4%, assuming a total luminosity of 300 fb~^ 

With this precision, systematic errors might become the limiting factor for the accurate 
determination of ^s- We have identified and estimated the following dominant sources of 
systematic uncertainties. 

The extraction of a cross section for a hard production process at a hadron collider 
depends on the knowledge of the parton distribution functions (PDFs) of the protons in 
the initial state. It is likely that our knowledge of PDFs for processes with large transverse 
momentum will be improved by data taken at the LHC itself. Nevertheless, here we adopt 
a conservative approach and estimate the uncertainty resulting from PDFs by comparing 
results for different CTEQ PDFs releases (CTEQ6M, CTEQ6D, CTEQ5M1, CTEQ5L and 
CTEQ6L1) [17] that are available today. It is found that this variation leads to an uncertainty 
of about 1% in the extracted Yukawa coupling. 
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Error source 


S[9s/9s] 


Method 


ABC 


LHC statistics for 300 fb^^ 


0.6% 3.0% 4.9% 


PDF uncertainty 
NNLO corrections 
Mass measurements Arriq = 10 GeV 
Am~g = 12 GeV 


1.4% 3.0% 3.0% 
2.0% 3.2% 3.2% 
2.0% 1.2% 1.2% 
1.1% 2.2% 2.2% 


Total 


3.4% 5.9% 7.3% 



Table 3: Combination of statistical and systematic errors for extracting the SUSY-QCD 
Yukawa coupling from like-sign dilepton events including jet transverse-momentum distribu- 
tions at the LHC. In Method A the complete information on squark branching ratios from a 
multi-TeV e~ collider is used for an SPSla-type scenario. Method B is based on LHC data 
potentially alone, but supported eventually by 1 TeV e~^e~ collider data, and combined with 
moderate assumptions on new heavy iso-scalar singlets. Method C allows for more complex 
scenarios in the neutralino sector than encoded in SPSla-type reference points, allowing 
again for the presence of new Standard Model singlets. 

Furthermore, the precise analysis of the squark and gluino cross sections relies on accu- 
rate squark and gluino mass measurements and on theoretical calculations of higher-order 
corrections. For the sparticle mass errors, we have adopted the results of Ref. [19], where an 
uncertainty of about 10 GeV is expected for the L-squark masses and about 12 GeV for the 
gluino masses. On the theoretical side, the next-to-leading order corrections are known [6] 
and have been included in the analysis. The uncertainty of the missing 0(as) contributions 
are estimated by varying the renormalisation scale of the 0{as) corrected cross section be- 
tween m/2 < Q < 2m, where m is the average mass of the produced sparticles (squark or 
gluino). 

The error estimate does not include the uncertainty for the squark branching ratios 
which are presumed to be small. In practice, due to the large squark masses in our reference 
scenario, they could either be measured at an e~^e~ linear collider with energy beyond 1 
TeV or, based alternatively on the closure of the MSSM, one could reconstruct the neu- 
tralino/chargino sector at a 500 GeV linear collider with high precision at the per-cent to 
per-mil level [3, 18] and derive the squark branching fractions from this information. 

The effect of the systematic errors and the combined error, obtained by adding all indi- 
vidual error sources in quadrature, are listed in Table [3] (Method A). It turns out that the 
total error, 6[gs/gs] = 3.4% is dominated by systematic effects. 

With a luminosity upgrade of the LHC, yielding a total luminosity of 1 ab~^, and as- 
suming that the mass and PDF uncertainties will improve with higher luminosity, this error 
could be reduced to 2.5%. 
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4 Towards model-independent analyses 



The procedure of determining the Yukawa couphng gs as outhned in the previous section 
can be expanded towards a more general analysis. Though the problem is too complex to 
allow for a completely model-independent analysis, methods can nevertheless be exploited 
which are based only on the pre-determined chargino/neutralino system, as expected to be 
known from LHC and/or e+e" collider measurements, combined with natural assumptions 
on the structure of theories beyond the minimal supersymmetric extension of the Standard 
Model MSSM. Central to this approach is the circumvention of external measurements of 
heavy squark-decay branching ratios which, for general expectations of the spectrum, can 
only be performed at a multi-TeV lepton collider. 

4.1 MSSM plus heavy iso-singlets 

First we will analyse a scenario in which the MSSM is expanded by a set of heavy iso- 
scalar particles, with little mixing between standard and new states, as could be expected 
in extended Higgs and gauge theories, see e.g. Refs. [20]. 

Proceeding from the same scenario as analysed in section [31 we observe that the charginos 
and neutralinos are light enough to be produced as diagonal and mixed pairs in e~^e~ collisions 
below 1 TeV, xtxj {hj = ^y'^) and XiX'j {hj = 1; • • • !4). Masses and mixing parameters 
can be measured in e~^e~ experiments, cf. Refs. [3], possibly even by observing only a 
subgroup of the lightest particles [21]. They determine the SU(2), U(l) gaugino and higgsino 
decomposition of the charginos and neutralinos. Thus the partial widths r(g — > x) can be 
predicted from these experimental results. Moreover, the decay branching ratios BR{x 

can be measured directly. In many scenarios even the observation of relatively large 
mass differences among the neutralinos and charginos [possible at LHC already] is sufficient 
to suggest strongly that the mixing in the neutralino and chargino sector is very small and 
thus can be neglected for the purposes here. This is the case, for instance, in the reference 
scenario studied in this paper. 

While the partial widths r{q —>■ x ^ i) can be completely determined in this scenario, 
the total widths r(g), and as a result the branching ratios, remain unknown since decay 
channels to heavy particles beyond MSSM, in particular to heavy iso-singlets, may be open. 
However, since R-squarks do not decay to leptons in scenarios where the lightest neutralino 
is predominantly bino, the only unknown quantities are the branching ratios BR{ul x^) 
and BR{dL — x~) to which all other branching fractions can be linked once the x mixing 
matrices are known; according to the Wigner-Eckart theorem: 

BR{uL^X2) = ^fuu{N^,N^)BR{uL^xt) 
BR{dL^xl) = \fM{N^,N^)BR{dL-^Xi) (10) 
The coefficients fuu,dd are ratios of the chargino A^='= and the neutralino mixing matrix 
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elements, with ^ being the ratio of the standard Clebsch-Gordan coefficients squared. For 
pure xf, X2 iso- vector winos, the relations simplify to fuu,dd 1- 

This could be the platform for a 3-parameter analysis, including the unknown branching 
ratios and the strong Yukawa coupling. However, the two branching ratios above are related 
if squark decays beyond xf and Xi involve only new iso-singlet particles with universal 
couplings to L-squarks, as may be expected, on quite natural grounds, outside the Standard 
Model. In the same notation as above, the relation reads 

BR{dL ^ XI) = fduiN^, iV°) BR{UL ^ xt) (H) 

[with fdu — i> 1 in the wino limit] . In this scenario the analysis of section [3] repeats itself 
except for leaving free the overall normalisation, governed by the unknown branching ratio 
BR{ul — > Xi)- In this case, the SUSY-QCD Yukawa coupling can only be extracted from 
ratios of the cross sections or from distributions with varying transverse momenta, but not 
from absolute values of cross sections. 

By analysing the measured rates for events with 2, 3 and 4 jets, one can extract a ratio 
of the cross-sections for ^lQ^l and qi^g production. This has also been exploited for the model 
analyses in Ref . [7] . By comparison with the theoretical calculations for these cross-sections, 
this ratio can then be interpreted in terms of the coupling ratio Qs/ Qs- 

The robustness and precision of the Yukawa coupling determination can be improved 
by looking not only at the number of identified jets, but also their transverse momentum 
distribution. 

The results of the simulations are presented in Figs. [31 As apparent from the figure, both 
the distributions as a function of the of the third and the fourth jet can discriminate 
between the two processes Eqs. ([2]) and ([3]). Squark-pair production tends to peak towards 
low values of pt,3 and pt,4, since the additional jets are generated in this process only by 
gluon radiation from the initial or final state. The qg Compton production and the gg 
gluino-pair production, on the other hand, typically lead to larger values for pt,3 and j9t,4- 
Therefore, the lower ends of the Pt,3 or pt,4 spectra depend more strongly on variations 
of than the region of higher transverse momenta. As a consequence, the shape of the 
transverse momentum distributions is sensitive to the ratio (js/gs, independent of the total 
normalisation of the cross-section. 

The expected experimental precision for the determination of the SUSY-QCD Yukawa 
coupling can be estimated from a binned x^ fit- Assuming 300 fb~^ of integrated luminosity, 
the statistical error for is 3.3%. This has to be combined with systematic errors as 

given in Table [3] (Method B), generating a final error of 5.9%. The systematic errors were 
estimated in the same way as explained in the previous section. 

Since the overall normalisation is left as a free parameter in this analysis, the result can 
also be exploited to determine the branching ratio of squark decays to charginos. A fit to 
Pt,3, for instance, can then simultaneously be interpreted as information about the Yukawa 
coupling and the squark branching ratio, see Fig. HI 
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Figure 3: Results of simulations for (a) the distribution with respect to the pi of the third 
jet, and (b) the distribution with respect to the px of the fourth jet (with pt,3 > 50 GeV). 
In both plots the contributions from the processes in Eqs. to as well as the SM 
background are shown separately, while for the sum the statistical errors for 300 fb~^ is also 
shown. 
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Figure 4: Independent results for the SUSY-QCD Yukawa coupling and the overall 
branching ratio of L-squarks into a single identified lepton / = e, /i (normalised to the MSSM 
value) from a fit to the p^s distribution in Fig. \^(a). The plot shows the la statistical error 
contours corresponding to the fit method B. The black bar indicates the result from the fit 
method A, assuming that the branching ratio has been pre-determined. 

4.2 Generalised picture 

In the second more general step we drop the assumption that X2 is nearly wino-like but 
allow, as before, for a set of new heavy singlets under the Standard Model group. 

In this scenario, the neutralino sector becomes more involved. Significant mixing between 
bino and wino components in the neutralino X25 for instance, can create a sizable branching 
ratio of R-squarks ur, into leptons [considering, as before, the decay X2 ~^ '^~^'^~Xi: with 
one tau decaying leptonically, the other hadronically, leads to one identified lepton in the 
final state]. Furthermore, mixing remnants with the extra iso-singlet components eventually 
may make the picture even more complicated. 

However, it should be noted that neutralinos always generate the same number of pos- 
itively and negatively charged leptons. On the other hand, chargino decay chains lead to 
leptons in the final state whose charge is correlated to the squarks from which the charginos 
originate. Therefore, this scenario can still be controlled by measuring separately the cross 
sections for i'^i'^ and i^£^ production. All neutralino decay modes of L- and R-squarks can 
be combined into one contribution, leading to equal rates of and i~ in the final state signa- 
ture. Then, two branching ratios remain undetermined, BR{ui — * Xi) ^"^^ BR{dL — ^ Xi )> 
that control the difference between the rates for and . 

This analysis relies mainly on the assumption that the chargino sector remains unmodified 
with respect to the MSSM when the system is expanded. [Only additional new charginos 
that couple differently to ul and di would significantly alter that picture.] In any case, the 
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Figure 5: Determination of the SUSY-QCD Yukawa coupling and the L-squark branching 
ratios into charginos from a fit to thep^s distributions for positive (l^t^) and negative (£~£~) 
lepton pairs in the final state. The plots show the la statistical error regions corresponding 
to the fit method C. 

chargino sector can be brought under control from measurements at a low-energy lepton 
collider. Furthermore, for our analysis of the signature of two leptons, two or more jets and 
missing energy, we will obtain a decent signal only for the mass hierarchy 3> ^ 
m^± S> 'm^Q^. For different mass spectra, the selection cuts would effectively kill the SUSY 
signal. It is likely that a similar analysis can be performed in such a case by suitably adjusting 
the cuts, but a thorough investigation of all possible situations is beyond the scope of this 
paper. 

We repeat the analysis of the previous section, now taking as measurable observables 
the [binned] spectra of the 3rd and 4th jet for and £~£~ separately. The free 

parameters in the fit are the Yukawa coupling and the branching fractions BR{ul — ^ Xi) 
and BR{dL — > Xi )> addition the overall normalisation and the contribution from R-squark 
processes, mainly — *• q^^q^^ . [On the methodological side, this case illustrates most clearly 
the technical potential of varying the relative weight of direct squark-pair production versus 
the Compton process in determining the Yukawa coupling.] Since the closure of the minimal 
super symmetric Standard Model, in the form defined above, is a natural assumption, the 
scenario is expected to be quite general and comprehensive. Weighing this qualification 
properly, the separate measurements come very close to a model-independent 

determination of the Yukawa coupling ^s- 

Assuming a luminosity of 300 fb~^, a binned fit gives a statistical error of 5.7% for 
gs, see Table |3] (Method C). The individual squark branching ratios, on the other hand, are 
determined relatively poorly, see Fig. [5l The origin of this uncertainty is the fact that both 
the pt spectra as well as the ratio of / signal rates discriminate a priori hardly 

between -ul and c^l decays. The distinction is slightly enhanced by the different shapes 
of u- and d-quark distributions in the proton, though the fit results for BR{ui — ^ Xi) 
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and BR{dL Xi ) remain strongly correlated. The correlation however does not affect the 
accuracy of the Yulcawa coupling significantly. With a relative contribution of 12% to ++ and 

20% to final states, respectively, X2 pairs contribute moderately to the signal events, 

but they lead to almost the same transverse jet momentum distribution as the chargino 
channels. This explains the relatively precise determination of Qs, even when the squark 
branching fractions are not pre-determined. 

The systematic errors in the analysis method C are equivalent to Method B. In total, 
the error on the ratio Qs/ Qs can be expected close to 7.3% in this largely model-independent 
approach. With a luminosity-upgraded LHC, yielding 1000 fb~^, this error could be reduced 
to about 5%. 

5 Conclusions 

A crucial consequence of supcrsymmetry is the identity of gauge couplings with the corre- 
sponding Yukawa couplings. Specifically in supersymmetric QCD, the quark-quark-gluon 
coupling and the quark-squark-gluino coupling are predicted to be equal [up to radiative 
corrections for SUSY scales different from the electroweak scale, described by renormalisa- 
tion group running]. In the present analysis we have investigated how this relation can be 
tested through measurements at the Large Hadron Collider (LHC). 

After identifying a clean signal of same-sign leptons, two or more jets and missing energy, 
that can be easily separated from Standard Model backgrounds, we have elaborated how this 
signal receives contributions from various squark and gluino production processes. Assuming 
a sufficiently large mass difference between gluino and squarks, it was further shown how 
gluino and squark production can be distinguished by analysing the transverse momentum 
distribution of additional jets. Since the gluino and squark production cross sections depend 
in different ways on the SUSY-QCD Yukawa couphng, this coupling can be extracted by 
measurements of the cross sections and jet distributions. 

We have established analysis strategies that allow the extraction of the SUSY-QCD 
Yukawa coupling with varying degrees of model assumptions. In the most restricted case, it 
has been assumed that the complete decay pattern of the squarks in the MSSM is known, 
possibly from a high-energy e^e" collider. In this case one can perform a measurement of the 
total cross sections for squark-pair, squark-gluino and gluino-pair production at the LHC. 

In a second step, we have extended the analysis to the situation when the squark branch- 
ing ratios are not known, allowing for a more general structure of the theory, with iso-singlet 
extensions of the MSSM, for instance. For relatively small mixing in the neutralino sector, 
several useful relations between different squark partial widths can be established based on 
gauge invariance, leaving only one branching ratio, ui to be determined internally 

within this LHC measurement. In case of large mixing in the neutralino sector, one ends 
up in general with three free parameters that nevertheless can be constrained from LHC 
data such that the impact on the measurement of the SUSY-QCD Yukawa coupling remains 
under good control. 

For a specific scenario we have demonstrated the apphcation of these methods. We have 
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generated Monte-Carlo samples of the squark and gluino processes as well as the relevant 
backgrounds. After applying cuts to reduce the background, we have extracted the SUSY- 
QCD Yukawa coupling from appropriate fits to the simulated data. Including the most 
relevant expected systematic errors, we have found that the SUSY-QCD Yukawa coupling 
can be determined at LHC with a precision of a few per-cent, i.e. 3.4% in the most restricted 
method, and 7.3% in the most general method. 

While these results are very encouraging, it should be pointed out that the experimental 
analysis and the expected quantitative precision will depend on the concrete supersymmet- 
ric scenario, i.e. the sparticle mass spectrum. Nevertheless, our results demonstrate that 
precision measurements of the underlying structure of supersymmetric theories are feasible 
at the LHC. It would be interesting to refine these studies in the future by including the 
impact of detector effects. 
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